A family of many-body models which are exactly solvable analytically 



I. Fuentes-SchulleU 

Instituto de Ciencias Nucleares, Universidad Nacional Autonoma de Mexico, 
A. Postal 70-543, Mexico D.F. 04510, Mexico and 
Perimeter Institute for Theoretical Physics, 31 Caroline St N, N2L 2Y5 Waterloo, Canada 

P. Barberis-Blostein 

Centro de Ciencias Fisicas, Universidad Nacional Autonoma de Mexico, Av. Universidad s/n, Cuernavaca, Mexico 

We present a family of many-body models which are exactly solvable analytically. The models 
are an extended n-body interaction Lipkin-Meshkov-Glick model which considers spin-flip terms 
which are associated with the interaction of an external classical field which coherently manipulates 
the state of the system in order to, for example, process quantum information. The models also 
describe a two-mode Bose-Einstein condensate with a Josephson-type interaction which includes 
n-particle elastic and inelastic collisions. One of the models corresponds to the canonical two-mode 
Bose-Einstein Hamitonian plus a term which we argue must be considered in the description of the 
two-mode condensate. Intriguingly, this extra term allows for an exact and analytical solution of the 
two-particle collision two-mode BEC problem. Our results open up an arena to study many-body 
system properties analytically. 



Many-body systems are of great relevance in most ar- 
eas of physics. In particular, there has been increasing 
interest in studying the properties of many-body systems 
and learning how to manipulate them in order to im- 
plement quantum information processing 1]. Ion traps, 
NMR systems, optical lattices, spin chains, and many 
others have been investigated for this purpose Q- Un- 
fortunately, interesting many-body systems are rarely ac- 
cessible to purely analytic analysis. Exact solutions exist 
mainly for one dimensional systems and the higher di- 
mensional cases have to be treated numerically 3J. Nu- 
merical calculations are in practice limited by the grow- 
ing degrees of freedom of the system. In nuclear physics, 
the Lipkin-Meshkov-Glick (LMG) model was intro- 
duced as a toy model to study many-body properties. 
It has been studied extensively because its integrability 
allows for numerical analysis [5j and approximate solu- 
tions using the algebraic Bethe ansatz This model 
has relevance in quantum optics since it is related to the 
two- mode Bose-Einstein condensate (BEC) Q. In the 
study of BECs, multicomponent condensates are of main 
interest. The lack of multimode BEC models which have 
analytical solutions has restricted the understanding of 
such mesoscopic systems whose most intriguing property 
is that they show collective quantum behavior. 

In this Letter we present a family of many-body models 
which are solvable analytically. The family is character- 
ize by the integer parameter n which indicates de max- 
imum number of particles which interact in the system; 
the n- model describes the n-body, (n-l)-body..., and 2- 
body interactions of N spin-1/2 particles which addition- 
ally interact with a classical coherent field. The 2-modcl 
correspond to an extended LMG model which considers 
single spin-flip terms produced by the interaction with 
the field and the effects of particle interactions during 
the spin-flip process. This extension allows for an exact 



analytical solution of the model. 

The n-model also describes a two-mode Bose-Einstein 
condensate with a Josephson-type interaction. The 
canonical Josephson Hamiltonian considers elastic 
two-particle collisions and has no analytical solution. So 
far, the canonical two- mode BEC model has not been 
quantitatively verified by experiment although it quali- 
tatively describes some of the observed effects . 

Our exactly soluble n-model corresponds to a more 
complete two-mode BEC model since it considers all the 
features of the previously mentioned model, but in addi- 
tion considers all n-particle elastic and inelastic collisions. 
Indeed, these multi-particle collisions must be included 
in an realistic description of the condensate. It has been 
extensively pointed out by experiment that inelastic col- 
lisions are present in the BEC and play an important 
role in some systems. Moreover, higher order collisions 
are relevant beyond the dilute regime of BECs [1E3- 
So far only a couple of models incorporate non-elastic 
collisions in multi-mode condensates . Many-particle 
collisions have not been addressed at all in the theoreti- 
cal models, although it is known that they are physically 
relevant especially in the coldest phase of the condensate 
where the particle density is high |lf| . 

Ironically, an effort is purposely made in the labora- 
tory to suppress many-particle collisions and inelastic 
processes, order to allow for comparison with the exist- 
ing theoretical models ||. Here we show that including 
these processes in the theoretical description, the model 
becomes exactly soluble. 

We show, using the 2-model, that the evolution of the 
relative population of the condensate presents collapse 
and revivals of Rabi oscillations. We calculate the ground 
state of the system, and show that under certain circum- 
stances, the ground state is in a multiple macroscopic 
superposition of coherent states. 
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The analytical solutions of this family of models will 
allow for a deeper understanding of many-body proper- 
ties. 

We introduce our models by considering the family 
of Hamiltonians Hq = Y^i=o A%J l z where J z is some 
representation of the SU(2) angular momentum opera- 
tor in the z direction and A, are real constants. Since 
J z \j, m) — m\j, m) with j and m integers or half integers 
with m = —j, — j + 1, j — 1, j, the eigenstates of the 
Hamiltonian are \j, m) with energy £ r " = J2i=o Aim 1 . By 
applying U = e*^ J *e y which is the most general rota- 
tion of J z in the SU(2) algebra with angles <f> and to 
Hq we construct the family of n-models, 

n 

H n = U^H^U = Ai(uU z Uf. (1) 

The exact and analytical solution of these Hamiltonians 
is of course simply U'\j, m) with energy £™ . The integer 
parameter n defines the n-model by considering up to n 
powers of J z in the ffj Hamiltonian. In terms of J± 
./.. ± i J yj the 2-model is 

H 2 = A 1 (cos6»J 2 +sin6»(e^J + + e- i<#, J_)) (2) 

+ A 2 (cos 2 0J 2 + sin 2 0(e«*4 + e~ 2 ^ J 2 _ 

+ J+J- + J-J+) 

+ cos 6 sin 9(J„(e*+J + + e"^J_) + h.c)). 

The Hamiltonian can be written in terms of two bosonic 
operators [a' , a] = [b', b] = 1 through the Schwinger rep- 
resentation which relates the bosonic opertaors to the 
angular momentum ones in the following way: J z = 
a^a — b^b, J+ = a^b and J_ = ab^ . Note that the comuta- 
tion relations of the SU(2) operators are indeed satisfied, 
and the total number operator N = n a + nj, = a + b^ b 
is related to the total angular momentum by J = N/2. 
Therefore, by choosing different representations of the 
SU(2) operators one can vary the total number of par- 
ticles N. The 2-model in the two-mode representation 
is 

H 2 = A + 5uj{a)a - b%) (3) 
+ Xie^a^b + e-^ab^+Ua^tfab 
+ A(e 2i *a t a t &6 + h.c.) 
+ ^((a t a t a&-& t a t a&)e i<#, + /i.a), 

with A Q = A 2 (cos 2 8N 2 + sin 2 ON), Slo = Aicosfl, 
U = A 2 (l - 3cos 2 6), A = Aism6, fi = 2 A 2 cos sin 
and A = A 2 sin 2 9. The states are related by \j, m) = 
\n a = j + m) ® | nft = j — m). Note that the particles are 
indistinguishable and the model only accounts for how 
many of them are in a given state. We will devote the 
rest of this paper to point out the physical relevance of 
the family of Hamitonians given by Eq. Q . 

The n-model describes the n-body, (n-l)-body,..., and 
2-body interactions of N — a^a + b^b spin- 1/2 particles 



(with n < N) in the presence of a classical coherent field. 
We will first analyze the 2-model. The first term in Eq. 
describes the free energy of a^a spin- 1/2 particles in the 
spin up state and b^b in the spin down state with fre- 
quency difference Soj . The interaction between two spins 
has strength U and corresponds to a dispersive process in 
which spins exchange their state while total spin is con- 
served. Additionally, we consider the interaction with an 
external classical field that produces one spin to flip state 
with coupling constant A. The classical field could be 
an effective field due to the presence of another system, 
other degrees of freedom of the system or, more interest- 
ing, to an external experimenter manipulating the state 
of the system using a laser. This last situation would be 
necessary for manipulation quantum information in the 
system. 

Due to the interaction of the field with the system, 
there is also a probability, parameterized by /i, of hav- 
ing two spins flip their state. Since the Hamiltonian has 
a second order character, i.e. it considerers products of 
two and four creation and annihilation operators, one 
must consistently consider all possible second order phys- 
ical processes. Thus, we include the two-particle spin-flip 
term (A) and the term that describes a single dispersive 
process (p) taking place due to particle interaction while 
the laser produces a single spin to flip. 

Consider an experimenter using a laser to change the 
state of a single spin. In an idealized situation the laser 
acts only on one spin and no interaction between spins is 
present. But in a realistic situation the field cannot be 
directed only to one spin and there is a probability that 
two spins change their state. One must also consider 
the possibility having a single dispersive process while 
the laser acts on a single spin due to spin interactions. 
Therefore the terms corresponding to fj, and A model a 
more realistic situation. It is remarkable that considering 
these extra terms allows for an analytical solution of the 
system. 

Now we can analyze what the n-model describes: all 
possible n-body, (n-l)-body,...,and 2-body interactions 
with n < N. Thus, Eq. QJ describes the possibility of n 
spin- 1/2 particles exchanging their state in such way that 
the total spin is conserved and considers a laser which 
causes m < n particles to flip state and all the possible 
dispersion terms accompanying this process. 

If no classical field is applied, the model simply consists 
of N spin-1/2 particles which interact in such way that 
the total spin is always conserved. The spin is conserved 
because no energy is provided to the system. 

Surprisingly, ignoring the single spin-flip term (A) and 
the term with a spin- flip plus single dispersion (/x), the 2- 
model corresponds to the LMG model of nuclear physics. 
The LMG model was constructed using products of two 
and four creation and annihilation operators with the 
purpose of creating a simple model for testing many- 
body properties. It has so far no physical realization 
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and it does not admit an exact analytical solution. Here 
we showed that considering an extension to the model, 
by considering consistently all possible products of two 
and four creation and annihilation operators, an exact 
solution is found. Our extension also provides the model 
with a consistent physical picture. 

Now lets focus on a closely related problem which does 
have a physical realization. The family of Hamiltonians 
Eq. JD is also of interest in BEC since it describes a two- 
mode BEC. Two mode BECs where first produced in 
JILAH3 and MIT0. The modes at, a and b\b corre- 
spond to atoms in two different hyperfine levels trapped 
by a magnetic potential with frequency difference 8u). 
The Hamiltonian considers the interaction with an ex- 
ternal laser to induce a Josephson-like coupling between 
the two modes with coupling constant A and phase <f> 
(l4|. This term corresponds to the annihilation of one 
particle in one mode and the creation of one particle in 
the other mode via the absorption or emission of a laser 
photon. The terms with four bosonic operators describe 
two-particle elastic and spin-exchange inelastic collisions. 
The elastic collisions have interaction strength U. The 
inelastic collisions have interaction strength fi when two 
particles of the same specie collide and one of them is 
transformed into the other via emission or absorption of 
a photon and interaction strength A when the collision 
transforms two particles of one type into the other. 

Note that by fixing 5uj, A and U the inelastic collision 
constants, A and /i, are determined. This is because in 
our model, the inelastic collisions are produced by the ef- 
fect of the laser on colliding particles. When the laser is 
applied to the condensate to induce transitions between 
the hyperfine levels, there is a probability of having a 
collision between two atoms and a laser photon. Since 
the laser is considered to be a classical field, this interac- 
tion gives rise to the two-particle inelastic collision terms. 
This physical relationship is mathematically expressed by 
the relationship between the coefficients. In our model 
the rate of elastic to inelastic collisions is given by 



A 



U 



A / A + 2Suj * 
2 \A\ -3<5cj 2 < 



(4) 



Ignoring the inelastic terms in Eq.Q, we find that our 
model coincides with the canonical Joscphson Hamilto- 
nian when the rate of collisions of same species is 
equal. The assumption of equal collision rates for the 
same specie is also made in J\ order to find approxi- 
mate and numerical solutions. Our model has the same 
number of free parameters as the canonical two-mode 
Hamiltonian. The only difference is that Hamiltonian 
in Eq.(|3J) includes inelastic collisions, which are usually 
present in real BECs 0. In magnetic traps, inelastic 
collisions are commonly suppressed to a large extent be- 
cause they give rise to atom losses. Spin exchange is 
the dominant loss mechanism in two-mode condensates 
[l5| . But in optical traps the particle loss due to spin 



exchange is negligible and there it is no longer necessary 
to suppress the process [l^. Including the correct rate 
of inelastic collisions in the Hamiltonian allows for an 
analytical solution which is simply U'\N,m). Note that 
U in the Schwinger representation is the two-mode dis- 
placement operator U = e^ b+ ^ ab ^ where £ = Oe 1 ^ . In 
|17j an author of this paper and collaborators proposed 
the Hamiltonian H2 to generate Berry phases in BECs 
but it was not understood then that this could indeed 
correspond to a solution of the two- mode BEC. 

We would like to emphasize that if for a specific phys- 
ical system, the relation in Eq.QJ between the rates of 
inelastic to elastic collisions does not hold or cannot be 
arranged by external manipulation, an analytical solu- 
tion cannot be found using our method for such a con- 
densate. Fortunately, in the laboratory the rate of elastic 
and inelastic collisions can be manipulated by applying a 
magnetic potential ,16] and it is possible to meet the ex- 
perimental values for the production ratios of those terms 
in a two-mode BEC. 

Let us now focus our attention on the 3-model Hami- 
tonian H3. This Hamitonian corresponds to all possi- 
ble three-particle and two-particle interactions including 
elastic and inelastic collisions between the same and dif- 
ferent specie. The n- model in Eq.JQ) describes a two- 
mode BEC where n-body, (n-l)-body,...,and 2-body in- 
teractions are considered. In the two-mode BEC n- 
particle collision terms are in principle present specially 
when the particle density is high. 

The canonical two-mode model 7] considers only two- 
particle elastic collisions and has no exact analytical so- 
lution. Commonly, the Bethe anzatz is used to find the 
ground and first excited state solution, or numerical work 
is needed. The model introduced here is more realistic, 
complete and has an exact analytical solution. It is pos- 
sible to analyze the whole spectrum and one need not 
restrict the attention only to the ground state. 

Due to the simplicity of our solution the ground state 
U'\N,mo) of H2 is trivially found by minimizing the 
energy — Aim + y^m 2 with respect to m. For 
A2 > 0, mo is the nearest integer to — Ai/(2A 2 ) or 
m = -AxN/lAt] when | -A 1 /(2A 2 )\ > N. For A 2 < 0, 
the minimum corresponds to mo = N if Ai < or 
mo = — N otherwise. The canonical two-mode BEC pre- 
dicts that the ground state of the condensate is, under 
certain conditions, a macroscopic superposition of two 
coherent states 0. In Fig. [T] we plot the relative popula- 
tion distribution for different ground states of the second 
order of our model and find that macroscopic superposi- 
tions can involve several coherent states for mo < 1000 
as shown Fig. 1(c) and Fig. 1(d) This difference must 
be due to inelastic collisions. 

The evolution of the relative population J z for a given 
initial state \ip(t = 0)) = Em=-JV C m W\N, m), with 
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FIG. 1: Ground state relative population distribution for 
2001 atoms. Different mo correspond to different intensities of 
the laser. Quantum superposition appears when mo < 1000. 



coefficients C m , is given by 

JV 

(jz) = -sin 6< ^ C m C m - 1 L r 



(5) 



-jV+l 



cos(0 + (£ ro _! - E m ) t)(N(N + 1) - m(m - 1)) 



The expectation value of J y which describes the evolution 
of the relative phase of the condensates is equal to J z plus 
a phase shift of ir/2. In Fig. [21 we plot the evolution of 
the relative population for the initial state \N, N) where 
the condensate consists of a single specie. The system 
presents Rabi-type oscillations with collapse and revivals. 
We are currently studying the effects of higher-order colli- 
sions in the oscillations. Interesting generalizations of our 



els replacing the SU(2) for SU(3) algebra. The first gen- 
eralization is preformed by applying a two-mode squeez- 
ing operator 5(a) = e aatb+_a ab to the Hamiltonian Hq 
which then has solution S\a)U^{0,(j})\j,m). The SU(3) 
model is found by applying the most general rotation in 
the SU(3) algebra to a polynomial in the SU(3) diago- 
nal generators. The bosonic representation of the SU(3) 
model describes a three mode (or spin-1) BEC. In prin- 
ciple the model can be extended to the 5(7 (n) algebra 
corresponding to spin-J condensates. A general method 
of finding models with exact analytical solution can be 
extrapolated from our model. For a given algebra find 
the diagonal generators, construct a polynomial in them 
and apply the most general rotation of the operators in 
the algebra to generate a new Hamiltonian. 

The understanding of many-body systems in dimen- 
sions higher than 1 has been limited by the lack of any 
realistic models with analytical solutions. The model we 
have introduced here allows for the first time an analyt- 
ical study of the n-body interactions of N spin-1/2 par- 
ticles in the presence of a classical coherent field and of 
a two-mode BEC with n-body elastic and inelastic col- 
lisions. The model extends the Lipkin-Meshkov-Glick 
model of nuclear physics and the canonical two-mode 
BEC models. The clear advantages of our model over 
these models include the possibility to study higher or- 
der interactions. Currently we study inelastic and many- 
body interactions in BECs and their effects in the phase 
transitions and entanglement properties of the system. 
We consider that the family of models that we have intro- 
duced opens an arena to study higher-dimensional many- 
body systems analytically. 

We would like to thank O. Castahos, J. Hirsh, O. 
Jimenez, J. Rogel-Salazar and P. Hess for interesting dis- 
cussions. 
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FIG. 2: Evolution of (j z ). The initial condition is \i>{Q)) = 
\N = 100, m = 100). Ai = 1, A 2 = 0.01, 6 = 3/2. 

family of models which are currently under study include 
a family of Hamiltonians with squeezing terms and mod- 
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